Landau-Khalatnikov Circuit model for Ferroelectric Hysteresis 



O 

o 

< 



I 

o 
o 



> 

00 

oo 

o 



X3 
O 

o 



X 



S. Sivasubramanian and A. Widom 
Physics Department, Northeastern University, Boston MA 02115 
(Dated: February 6, 2008) 

We present the circuit equivalent of the Landau-Khalatnikov dynamical ferroelectric model. The 
differential equation for hysteretic behavior is subject to numerical computer simulations. The size 
and shape of the simulated hysteretic loops depends strongly on the frequency and the amplitude of 
the driving electric field. This dependence makes the experimental extraction of the coercive electric 
field difficult. The bifurcation of the driven Landau-Khalatnikov model is explained in detail. 

PACS numbers: 77.80.Dj, 77.80.Fm, 77.22.Gm 



I. INTRODUCTION 

Many of the interesting properties of ferroelectric ma- 
terials are probed experimentally by measuring the po- 
larization response to a time varying electric field. These 
properties form the basis for important potential com- 
puter engineering applications |jl|, ||, |[ Q such as ferro- 
electric random access memory or ferroelectric memory 
field effect transistors. A common method to observe fer- 
roelectric polarization features is via the hysteresis loop. 
A sinusoidal electric field is applied to the ferroelectric 
sample while the polarization P and/or polarization cur- 
rent density J — (dP /dt) is continuously monitored. The 
observed data is then in principle utilized to determine 
the thermal remnant polarization Ps and coercive field 

Some of the partially understood experimental results 
are as follows: (i) No unique coercive field is evident 
during the hysteretic process of switching, (ii) When 
an alternating electric field is applied, a hysteretic loop 
is observed only within a limited frequency band and 
a limited range of amplitudes, (iii) Multiple hysteresis 
loops are commonly observed. 

Various approaches have been adopted to explain these 
strange characteristic features of the ferroelectric ma- 
terials. One approach is to create a equivalent cir- 
cuit models I, ll 0, I, I, |ll|, [l2[ |l3| that simulate 
the experimental data. The earliest among these mod- 
els are the Sawyer- Tower circuit and modifications 
thereof||l^, |l^ which are used to study the hysteretic 
properties of ferroelectric capacitors. Currently, many 
of the circuit models are computed using the SPICE 1 17 
simulations. Previous circuit models depend strongly 
upon which particular applications one is interested in 
analyzing. No single standard circuit has been put for- 
ward to explain the great variety of experimental re- 
sults available for different regimes of electric field and 
frequency. Our goal is to explore many of the exper- 
imentally observed properties employing a single cir- 
cuit model. The model is equivalent to the Landau- 
Khalatnikov dynamical equation for the ferroelectric po- 
larization. This non-linear circuit model produces a va- 
riety of previously unexpected hysteretic behaviors. 

In Sec. II we present the circuit equivalent of the 



Landau-Khalatnikov|18| |l9| , pO| , gl[ dynamical ferroelec- 
tric model. In Sec. Ill the differential equations which 
must be solved are written in the "dimensionless" form 
suitable for numerical simulations. In Sec. IV the depen- 
dence of the size and shape of the simulated hysteretic 
loops on the frequency and amplitude of the driving elec- 
tric field will be exhibited. Sec. V deals with the bifur- 
cation of hysteretic curves. The dynamical properties 
of the model contains a regime of broken symmetry and 
a regime of unbroken symmetry. The conventional hys- 
teric loops are in the symmetric dynamical regime. In 
the concluding Sec. VI future experimental prospects for 
the model are explored. 



II. THE LANDAU-KHALATNIKOV CIRCUIT 

The thermodynamic equations of state for a ferroelec- 
tric model are described by the energy per unit volume 
U{P,S) as a function of polarization and entropy per 
unit volume obeying 



dU = E. 



thermal 



■ dP + TdS. 



(1) 



In addition to the thermal electric field, a ferroelectric 
sample exhibits a dissipative electric field derived from 
an Ohm's law resistivity p. 
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The total electric field E = Ethermai + 'Eohm gives rise to 
the Landau-Khalatnikov dynamical equation of motion 



E 



dPJs 



dP 

'dt 



(3) 



Finally, the Maxwell displacement field D within a fer- 
roelectric material is given by 



D = eoE + P. 



(4) 



If the ferroelectric material is placed inside of a capac- 
itor, then the total charge on one capacitor electrode is 
given by the surface integral 



Qtot = (P D • dS. 

electrode 



(5) 
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Were the capacitor electrodes embedded in the vacuum, 
then the capacitor charge would be 

Qvac = eo / E • dS. (6) 

electrode 

In general, Eqs.(4), (5) and (6) imply 

Qtot = Qvac + Q (7) 

where 

Q= (f P • dE. (8) 

J electrode 

Eq.(7) describes the charges on two capacitors in paral- 
lel with the geometrical capacitance Co relates the charge 
to voltage ratio of the vacuum 

Qvac = CoV. (9) 

Eqs.(3) and (8) describe the Landau-Khalatnikov dynam- 
ics for the non-linear capacitor in FIG. 1 with voltage 
V(Q) connected in series with an Ohm's law resistor R; 
i.e. 



Co 



li h({Q,S) denotes the energy of the non-linear capaci- 
tance, then 

dU = TdS + VdQ. (11) 
The circuit version of Eq.(3) is that 

-(iX-(f)^ <'^) 

The circuit Eq.(12) is pictured in FIG. I. In the Landau- 
Khalatnikov circuit, the upper capacitor carries a charge 
Qvac = CqV while the lower non- linear capacitor carries 
a charge Q and voltage V determined by the thermody- 
namic Eq.(ll). The resistance R describes the dissipa- 
tion present when the polarization varies with time. 



III. PERIODIC VOLTAGE SOURCES 

Hyst(^rc^tic cycles are measured applying a time varying 
voltage of the form 

V{t) = Vo cos{ujt). (13) 
For the case of the Landau energy shown in FIG. 2, 




where the saturation polarization Pg determines 

Qs^<f Ps-d^, (15) 

J electrode 



V Q 




FIG. 1: The equivalent circuit to the Landau-Khalatnikov 
dynamical equation is shown in the above figure. The lin- 
ear capacitor is on the upper branch of the circuit and the 
non-linear capacitor is on the lower branch. The polarization 
energy U is stored in the non-linear capacitor. The ohmic 
resistor R describes the dissipation due to time variations in 
the ferroelectric polarization. The geometric capacitance Co 
is defined in Eq.(9). V is the applied voltage across the par- 
allel circuit. 
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FIG. 2: Shown is the Landau energy for a ferroelectric model. 
The energy stored in the non-linear capacitor is given hy U = 
(Q?/8C){1 - {Q/Qsr)y . 
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Eqs.(12), (13) and (14) read 




= Vocos{Lut). (16) 



In order to solve Eq.(16) numerically, let us introduce 
the dimensionless quantities 

e = ut, y = {Q/Q,), rj={2toRC)-^ 

and 

z={Vo/RioQ,). (17) 
Eqs.(16) and (17) now read 

+w(y'-i) = ^cos^ (18) 

One seeks solutions to Eq.(18) which are periodic 

y{9 + 2TT;z,r]) = y{0;z,r]). (19) 

For example, let us suppose that at an initial time zero 

y{0 = O;z,ri) = x. (20) 

After a numerical integration through one period of mo- 
tion, one then finds that 

y{e = 2Tr;z,r]) = g{x;z,r]). (21) 

Eqs.(18), (20) and (21) define the function g{x; z, if). It is 
possible to analytically compute Q{x; z, rf) in two limits: 

Ymig{x;z,ri) =x, (22) 

T7— »0 



unique periodic solution for y{9; z,i]). A hystcrctic loop 
located in the {Q, V) plane may be found by eliminating 
the time parameter t in a one period interval < t < 
{2Tr/w) from the parametric equations 

Q = Q^y{e = ojt; z, ri) and V = RjjQsZ cos{ojt). (25) 

In FIG. 3, we have plotted two hysteretic loops corre- 
sponding to the same frequency but with different driv- 
ing voltage amplitudes. The inner smaller hystcrctic loop 
corresponds to a smaller driving voltage than that of the 
outer larger hysteretic loop. In FIG. 4 we have plot- 
ted two hysteretic loops corresponding to the same volt- 
age amplitude but with different frequencies. The inner 
smaller hysteretic loop corresponds to a lower frequency 
than that of the outer larger hysteretic loop. 

Since the large amplitude and high frequency loops 
have a larger enclosed area than do the small amplitude 
and low frequency loops, it is not at once evident how to 
extract the coercive voltage directly from hysteretic loop 
data. This is a well known problem in attempting to 
measure the coercive electric field forcing a polarization 
flip. The nature of the problem becomes evident from 
a study of the effective circuit in FIG. 1. The theoreti- 
cal "coercive voltage" Vc is that voltage across the lower 
branch non-linear capacitor present at the time when the 
polarization flips. The experimental "coercive voltage" 
Vc ''^^^ = RIc + Vc includes the voltage across the resis- 
tor. It is not an easy matter to separate the dissipative 
voltage from the thermal capacitor voltage. 

In order to understand the "area" enclosed by the hys- 
teresis loop in the (Q, V) plane, one notes that the ther- 

2 I 1 1 



lin.g(x;.,.)-f . (23) 
V V 1 + a;^(exp(47r?7) - 1)/ 

A sufficient condition for existence of periodic Eq.(19) 
solutions of Eq.(18) is that a; be a fixed point (solution) 
of the equation 

x = g{x;z,vi). (24) 

If Eq.(24) has a unique solution for x, than there exists a 
unique periodic solution of Eq.(18). If Eq.(24) had more 
than one solution for x, then there will be (in general) 
more than one periodic solution to Eq.(18). The number 
of solutions for x depends on the values of parameters in 
the {z,ri) plane. 

IV. HYSTERESIS CURVES 

In this section we consider the region in the {z, rj) plane 

for which the following conditions hold true: (i) Eq.(24) 
has a unique solution for x and (ii) there there exists a 




_2 I 1 1 

-1.5 1.5 

(CV/QJ 

FIG. 3: Shown are two different hysteretic loops simulated 
for identical frequencies but differing amplitudes of the ap- 
plied voltage. The smaller inner hystcrctic loop corresponds 
to z = 20.0 and r] = 25.0. The larger outer hystcrctic loop 
corresponds to z = 200.0 and r? = 25.0. Large applied ampli- 
tudes imply high apparent values for the experimental coer- 
cive voltages. 
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modynamic work performed by the non-linear capacitor 
during one cycle is given by 

-MK (cycle) = / VdQ, (26) 

J loop 

which is also (via the first law of thermodynamics) the 
heat dissipated by the resistor R during the cycle. The 
area enclosed by the hysteretic loop is large or small, 
respectively, when the dissipated heat in the resistor is 
large or small. 

In order to understand the role of dissipative heating 
in more detail, multiply Eq.(lO) by the current 

-(f) <-) 

and integrate the result over time for one cycle; i.e. 
/ VIdt = I (^)dt+(f RlHt. (28) 

J loop J loop \dQj \ dt J Jlggp 

The first term on the right hand side of Eq.(28) obeys 
(f dU = 0. (29) 

J loop 

We then see directly that the work done by the non-linear 
capacitor during one cycle must be dissipated as heat in 
the resistor; i.e. Eqs.(26), (27), (28) and (29) imply 

-VF(cycle) = (f VIdt = (f Rl'^dt. (30) 

J loop J loop 

Thus, a large applied voltage amplitude and/or a large 
applied voltage frequency yield a large area loop because 
more heat is dissipated per cycle by the resistor. 

2 I 1 1 
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FIG. 4: Shown are two different hysteretic loops simulated 
for identical amplitudes but differing frequencies of the ap- 
plied voltage. The smaller inner hysteretic loop corresponds 
to z = 100.0 and rj — 80.0. The larger outer hysteretic loop 
corresponds to z = 3.0 and tj = 2.25. Large applied frequen- 
cies imply high apparent values for the experimental coercive 
voltages. 



V. DYNAMICAL BROKEN SYMMETRY 

From the parity symmetry of the energy equation of 

state U{—Q,S) = 14{Q,S), the are two possible thermal 
equilibrium minimum energy states at Q± = ±Qs in the 
absence of a driving voltage. For a small driving ampli- 
tude Vq in Eq.(13), we then expect two possible periodic 
charge response functions Q± {t + [2n juS) = Q,±(t)\ 

Q^(t)^Q,_(t) =4> Broken Symmetry . (31) 

On the other hand, for a large driving amplitude Vq in 
Eq.(13), we expect a unique periodic charge response 
functions Q{t^ (27r/a;) = Q(t); 

Q_^(t) = Q_(t) =^ Restored Symmetry . (32) 

In terms of the dimensionless variables in Eq.(17), 

z = VqI{RujQ.s) and r]=\l<(lu)RC), (33) 

one expects the two regimes to be described by a dynam- 
ical bifurcation fimction B{r\) 

z < B{ri) =^ Broken Symmetry , (34) 

z > B{r}) Restored Symmetry . (35) 

We have numerically computed B(Tf) and the results are 
plotted in FIG. 5. The regions in the {z, rf) plane are ex- 
hibited corresponding to the symmetric and broken sym- 
metry dynamical phases. In the broken symmetry phase, 
there are two possible hysteretic curves; Q+{t) is local- 
ized in the neighborhood of +Qs and Q-{t) is localized 
in the neighborhood of — Qg. 
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FIG. 5: Shown is the phase plane for the Landau-Khalatnikov 
model with the dynamical bifurcation curve z = B{r]). If 
z < B{r]), then parity symmetry is broken. li z > B{ri), then 
the parity symmetry is restored. 
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For a more clear picture of the bifurcation |2j] 
boundary curve at z — B{ri) we show a sequence of hys- 
teretic loops in FIG. 6. For r/ — 0.10, we have chosen 
four values for the amplitudes z. Explicitly 

< Za = 0.1 < Zb = 0.7 < B{t] = 0.1) 

for illustrations of broken symmetry, and 

B{t] = 0.1) < Zc = 1.0 < Zd = 20.0 

for illustrations of unbroken symmetry. The two ovals in 
both FIG. 6(a) and FIG. 6(b) correspond, respectively, to 
Q±{t) for the values Za and Zb. The symmetry is broken 
because only one of the ovals (Q+ or Q_ but not both) 
will appear in an experiment described by the model. 
The oval in the unbroken symmetry phase corresponding 
to Zc is unique but does not appear very similar to con- 
ventional hysteretic curves. By increasing the amplitude 
to Z(j > Zc a conventional looking hysteretic loop appears 
in FIG. 6(d). 

VI. CONCLUSIONS 

It has been shown that the Landau-Khalatnikov model 
of ferroelectric hysteresis has a very rich dynamical struc- 
ture. Only a small part of the symmetric dynamical 
phase had been reported in standard literature. The 
conventional looking hysteretic curves are present the- 
oretically and experimentally only in the regime of high 



driving amplitudes. The region in the neighborhood of 
the bifurcation into the broken symmetry phase has re- 
ceived little or no experimental attention. The numerical 
simulations of this work indicate what is to be expected. 
Experimental measurements of the dynamical bifurcation 
curve z = B(r]) of ferroelectrics would be of great inter- 
est in checking the validity of the Landau-Khalatnikov 
model. 

Further insights on how to measure the coercive volt- 
age can be obtained by driving the circuit in FIG. 1 with 
a voltage source which has both DC and AC components 

Vit) = Veo.t + ViCosiu;t). (36) 

The DC voltage component Vext can be put into the en- 
ergy function via 

U{Q, S) ^ U{Q, S; Vext) = W(Q, S) - V^xtQ, (37) 
leaving Eq.(12) in the form 

V, cos{U) ^ (38) 

The advantage of employing a DC component is now 
evident. The symmetry breaking of the energy in FIG. 2 
will now be controlled by varying the DC voltage. The 
DC offset will allow for precise measurements of the co- 
ercive field. 
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FIG. 6: Shown are four loop plots choosing 77 = 0.1. (a) 
For the small amplitude B{ri = 0.1) > Za = 0.1, there are 
two possible small oval loops. Only one will be realized in a 
particular experiment, (b) As the voltage amplitude is in- 
creased to B(r] — 0.1) > Zb = 0.7 > Za the ovals grow 
large but are still representative of two different loops, (c) 
In the symmetric phase Zc — 1.0 > B{rj — 0.1), only a single 
loop is allowed which does not yet look very close to con- 
ventional hysteretic loops, (d) For large voltage amplitude 
Zd ~ 20.0 > Zc > Bijj = 0.1), the conventional hysteretic 
loop is recovered for the model. 
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